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ABSTRACT 

This paper deals with the approximation the solution of singular Riccati differential equations using the reproducing kernel Hilbert space scheme. The exact 
solution u(r) is represented in the form of series in the spacew_22 [0,1]. In the mean time, the n-term approximate solution u(r) obtained and is proved to 
converge to the exact solution u(r). Some numerical examples have also been studied to demonstrate the accuracy of the present method. Numerical experiments 
are performed to confirm our theoretic findings. 
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1. INTRODUCTION: 
In this work, we consider the following quadratic Riccati differential equation with singularity in reproducing kernel space 


PW)"(Y) = Q(*NUW7(r) + Ru) + S(r), uO) = & 0<r<il, (1) 
where P(r), Q(r), R(r) and S(v) coefficients are continuous real functions, perhaps P(0) = 0 or P(1) = Oand u(r) € w?[0,1] 


Many physical phenomena such as theory of thermal explosions, studies of electro-hydrodynamics, nuclear physics, gas dynamics, 
chemical reaction, studies of atomic structures and atomic calculations can be modeled by singular differential equations. Due to 
this reason, in the literature there exist numerous methods dedicated to the singular problems, for example see [1—3]. 


The theory of reproducing kernels, firstly has been presented based on the S. Zaremba works in [4] for solving boundary value 
problems with harmonic and bi-harmonic functions. In recent years, reproducing kernel methods emerged as a powerful framework 
in numerical analysis, differential, integral equations, probability and statistics [5—7]. 


Furthermore, extensive works have been proposed and discussed based on their theory for the numerical solutions of singular 
problems. 


For example, using operator transformation considered the solution of Cauchy singular integral equation by reproducing kernel 
Hilbert space (RKHS) scheme in [8]. Also, in [9] has been used from RKHS method to approximating the solution of the second 
kind weakly singular Volterra integral equation on graded mesh. 


In [10-17] have been solved various classes of singularly boundary value problems. In [18—20] the authors has been implemented 
the RKHS method for solving singularly perturbed differential-difference equations. 


A brief outline of this paper is as follows: 

In section 2, RKHS is introduced. In section 3, the analytical solution for (1) in the space w[0,1] is introduced. Also, an iterative 
method to solve (1) numerically in the space w?[0,1] is described. In section 4 the numerical examples are discussed. 

Then we end with a brief conclusion in section 5. 


2. REPRODUCING KERNEL SPACES: 
In order to solution of (1), we construct a reproducing kernel functions. 


Definition 2.1. ((21]) Let E be a nonempty abstract set and C be the set of complex numbers. A function K: E x E— Cisa 
reproducing kernel of the Hilbert space H if 


e for eacht € E, K(,t) € H, 
e For eacht € Eand w € H, (w(.), K (, t)) = w(t). 
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In Definition 2.1, second condition called the reproducing property, the value of the function yw at the point ¢ is reproducing by the 
inner product of y(t) with K(.,f). 


Remark 2.1. A Hilbert space H of functions on a set E is called a RKHS if there exists a reproducing kernel K of H. That is, a 
Hilbert space which possesses a reproducing kernel is called the RKHS. 


Definition 2.2. The inner product space w3 [0,1] is defined as w3[0,1] = {x(t)|x(t),x‘(t), is absolutely continuous, x(0) = a and x"(t) 
€ L2[0,1]}. 


The inner product w3[0,1] is defined by 
1 


< x(t), y(t) > w[0,1] = x(0)y(0) + x(1)y(1) + | x (t)y (t)dt, (2) 


0 
and the norm |\|x||w2[0,1] is denoted by ||x||w2[0,1] = /<x,x > w2[0,1], where x,y € w2[0,1] 
Theorem 2.1. ( [22]) The space w3 [0,1] is a reproducing kernel space. 

That is, for any x(t) € w3[0,1] and each fixed t € [0,1] there exists R,(t) € 

w2[0,1] , such that < x(t),R,(t) > w3[0,1] = x(z) . The reproducing kernel R,(t) can be denoted by 


- ((z — 1)t? + z(z? — 32 + 8)),t <z, 


R,(t) = 2 (22(¢ - 1) + t(t? - 3t+8)),t > Z. 


(3) 


Definition 2.3. The inner product space w2[0,1] is defined as w2[0,1] = {x(O|x(0), is absolutely continuous, x(0) = a and x'(t) € 
L7[0,1]}. The inner product w2 [0,1] is defined by 


< x (t), y(t) > w3[0,1] = x(0)y(0) + J x’ Oy’ (Oat, (4) 


and the norm ||x||w3[0,1] is denoted by||x|]w3[0,1] = /<x,x >w3[0,1] , where 
xy € w3[0,1]. 


In [5] proved that w3[0,1] is a complete reproducing kernel space and its reproducing kernel is 


= _1+t+t, t<zZ, 
R(t) = ens t>z, () 


3. THE ANALYTICAL SOLUTION OF (1) 
In this section, the solution of (1) is given in the reproducing kernel space w3[0,1] 


In (1), it is clear that £:w3[0,1] — w4[0,1] is a bounded linear operator. 


Put y.(z) = Rzs(z) and ¢,(z) = L*y;(z), where L* is the adjoint operator of L. The orthonormal system {9, (z)} 
[ee) 
s=1 


[oe] 


ee: w3[0,1] can 


be derived from GramSchmidt orthogonalization process of {9, (z)} 


@,(z) = > Xsi Gi (z), (“> 0,5 = 1,2,...) 
i=1 


[oe] 


Theorem 3.1. ([22])Let {zs}, _ 4 


is dense on [0,1], then {0, (2)}. iS 1 is the complete system of w$[0,1] and $,(z) = LtRz(t)|t=zs. 


Theorem 3.2. If {Zs}. is 1 is dense on [0,1] and the solution of (1) is unique, then the solution of (1) satisfies the form 


v(z) = Yih, & sif (zi, v(z))O.(z) (6) 





Research Paper E-ISSN No : 2454-9916 | Volume: 6 | Issue: 10 | Oct 2020 





Proof. From Theorem 3.1, one can easily see that {, (z)}. ie 1 is the complete orthonormal basis of w#[0,1]. Remark that from 


reproducing property (v(z),7s(z)) = v(zs) for each v(z) € w3[0,1], so we have 
v(2) =) < v@),B.(@) > wH[0,1]8,(2) 
s=1 


Ss 


= ey « si < v(z),L£5,(z) > w3[0,1]8,(z) 


i=1 
foe) s 


= 3 >» « si < Lv(z),x5(z) > w3[0,1]@,(z) 


s=1 i=1 
o s 


: > x si < f(z,,v(z)), x5(z) > w2[0,1]0,(z) 


s=1 i=1 


= . y « sif (z;, v(z;)),9,(2), 


s=1i=1 
and the proof of the theorem is complete. 


3.3 The implementation scheme: 
Equation (6) can be denoted by 


v2) = ) AG.) 
s=1 
Where A, = ij-1 “si f (Z;, v(Z;). Let zi = 0, it follows that f(zi,0(z1)) is 


known. Considering the numerical computation, we put vo(z1) = o(z1) and define the n-term approximation to v(z) by 


UplZ) = Fie BLO.(Z) (7) 


where the coefficients Bs and ¢,(z) are given as 
By = %y fa vo(Z1)), 
v,(z) = B,9,(2), 


2 


B= > Xi f (Zi, vi-1@)) 


v,(Z) = eek B,O,(z), 
Bid = Vien-1 %2i f (zi, vi-1 (Z)) 
vp (Zz) = Yo1 B,G,(z). 


In the iterative process of (7), we can guarantee that the approximation 0,(z) satisfies the initial condition u(O) = ¢. On the other 


hand, the approximate solutionUn (z) can be obtained by taking finitely terms in the series of v,(z) as follow, 


vn (z) = D1 Dies Xsi f Zi Vn-1(Z:))O5(Z) (8) 


4. NUMERICAL SIMULATIONS 
In this section, the scheme in the paper will be applied to four numerical examples. All of the computations have been performed 
by using the Matlab R2010a. Results of each example are compared with exact solution. 


Problem 4.1. In equation (1), if P(r) =r, Q(r) = 1, RW) = —Vr and S(”) =~-] then true solution is 
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RKHS method, taking with the reproducing kernel function Rz(t) on [0,1], the numerical solution uy, Ncr) computed by (8). The 
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2vr —3 
= 
—2,+4vr-3 
s-1 
u =1,, = ——,s =1,2,..., 
0 Lr No 1,2 N 


numerical results at some selected gird points for N = 51 and n= 7 are given in Table 1. 


Table 1: Numerical results for problem 4.1. 





r 
0.1 
0.2 


True solution 
1.223480959455358 
1.306879269928204 
1.351601504414804 
1.366020440416440 
1.353553390593274 
1.316983583242455 
1.259474520209441 
1.184736379760737 
1.096856471680698 
1.000000000000000 


Approximate solution 


1,223480953436179 
1.306879257935626 
1.351601534632646 
1.366020491746745 
1.353553470343467 
1.316983693532165 
1,259474421125751 
1.184736316453742 
1.096856438463759 
1.0000000253 19537 


Problem 4.2. Consider the following singular equation 


Absolute error 


6.0192 x 10-2 
1.1993 x 1075 
3.0218 x 107° 
5.1330 x 107° 
7.9750 x 10-5 
1.1029 x 1077 
9.9084 x 107° 
6.3307 x 107-* 
3.3217 x 1075 
2.5320 x 1075 


Relative « error 
4.9197 x 10-9 
9.1765 x 107° 
2.2357 x 107° 
3.7577 x 107° 
5.8919 x 10-5 
8.3744 x 107° 
7.8671 x 1075 
5.3436 x 1075 
3.0284 x 1075 
2.5 


5320 x 107° 





(—r)w(r) = w(r)tu(r), ~—u(0) = -5, 0 <r <1, with true solution 
u(r) = —5 
s-1 
t= Vu ri = 1,2,...,N 





Figure 1: Comparisons of approximate solution with the exact solution (left) and the absolute errors with the relative errors of 
problem 4.2 (right). 
Problem 4.3. Consider the following singular equation 
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u(r 
u(r) — vru?(r) + uae u(0)=1, O<r<il 
Vr 
() : 
LL) > jj —]}. 
1+2(r—vr)+ e-2vr 
by RKHS method, taking r, = —,s = 1,2,...,N on [0,1]. The numerical results at some selected gird points for N= 53 andn= 
9 are given in Table 2. 
Table 2: Numerical results for problem 4.3. 
r True solution Approximate solution Absolute error Relative error 


0.1 1.820117633910813 1.820117625479621 8.4312 x 10-9 4.6322 x 10-° 
0.2 2.187191844632269 2.187191765623859 7.9008 x 10-° 3.6123 x 107° 
0.3 2.383944752470398 2.3839447257395 16 Pa ye eh bere 


0.4 2.446922176511700 2.446921486935627 6.8958 x 1077 2.8181 x 1077 
0.5 2.412827055573784 2.412826748392754 3.0718 x 10-7 1.2731 x 1077 
0.6 2.316890586034218  2.316890058361949 5.2767 107" 2.2775 x 1077 
0.7 2.187458502978635  2.187457726387531 7.7659 x 107" 3.5502 x 1077 
0.8 2.044368584151155 2.044368509754165 7.4397 x 10-° 3.6391 x 107° 
0.9  1.900063469836569 1.900063425736172 4.4100 x 10-5 2.3210 x 107° 
1.0 1.761594155955765  —«:1.761594138936148 1.7020 x 10-° 9.6615 x 107° 
Problem 4.4. Consider the following singular equation 
(1—et)u'(r) = u(r), u(0) = = 0<r<1, with true solution 
1 
u(r) er+i-e 
Using our method, take r, = —,s = 1,2,..,N N = 35. andn=7. The numerical results are given in Figure 2. 





0 0.2 04 0.6 08 1 
t 
Figure 2: Comparisons of approximate solution with the exact solution (left) and the absolute errors with the relative 
errors of problem 4.4 (right). 





5. CONCLUSIONS: 
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In this paper, the reproducing kernel Hilbert space scheme was implemented for solving a singular Riccati differential equations. 
This confirms the validity of the present method and it is efficient, accurate and reliable for singular Riccati differential equations. 
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